New Pivot Selection for Sparse Symmetric Indefinite Factorization
نویسندگان
چکیده
We propose a new pivot selection technique for symmetric indefinite factorization of sparse matrices. Such factorization should maintain both sparsity and numerical stability of the factors, both of which depend solely on the choices of the pivots. Our method is based on the minimum degree algorithm and also considers the stability of the factors at the same time. Our experiments show that our method produces factors that are sparser than the factors computed by MA57 [9] and are stable.
منابع مشابه
A study of pivoting strategies for tough sparse indefinite systems
The performance of a sparse direct solver is dependent upon the pivot sequence that is chosen during the analyse phase. In the case of symmetric indefinite systems, it may be necessary to modify this sequence during the factorization to ensure numerical stability. Delaying pivots can have serious consequences in terms of time as well as the memory and flops required for the factorization and su...
متن کاملA New Pivot Selection Algorithm for Symmetric Indefinite Factorization Arising in Quadratic Programming with Block Constraint Matrices
Quadraticprogrammingis a class of constrained optimizationproblem with quadratic objective functions and linear constraints. It has applications in many areas and is also used to solve nonlinear optimizationproblems. This article focuses on the equality constrained quadratic programs whose constraintmatrices are block diagonal. Using the direct solution method, we propose a new pivot selection ...
متن کاملCompressed threshold pivoting for sparse symmetric indefinite systems
A key technique for controlling numerical stability in sparse direct solvers is threshold partial pivoting. When selecting a pivot, the entire candidate pivot column below the diagonal must be up-to-date and must be scanned. If the factorization is parallelized across a large number of cores, communication latencies can be the dominant computational cost. In this paper, we propose two alternati...
متن کاملTowards a stable static pivoting strategy for the sequential and parallel solution of sparse symmetric indefinite systems
We consider the direct solution of sparse symmetric indefinite matrices. We develop new pivoting strategies that combine numerical and static pivoting. Furthermore, we propose original approaches that are designed for parallel distributed factorization. We show that our pivoting strategies are numerically robust and that the factorization is significantly faster because of this static/numerical...
متن کاملImproving the stability and robustness of incomplete symmetric indefinite factorization preconditioners
Sparse symmetric indefinite linear systems of equations arise in numerous practical applications. In many situations, an iterative method is the method of choice but a preconditioner is normally required for it to be effective. In this paper, the focus is on a class of incomplete factorization algorithms that can be used to compute preconditioners for symmetric indefinite systems. A limited mem...
متن کاملذخیره در منابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
عنوان ژورنال:
- CoRR
دوره abs/1601.06812 شماره
صفحات -
تاریخ انتشار 2016